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Results from a study of the Nambu-Jona-Lasinio model on the lattice 
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The main results of our analysis of the two flavor Nambu-Jona-Lasinio model with SU(2) x SU(2) chiral 
symmetry on the four-dimensional hypercubic lattice with naive and Wilson fermions are presented. Large 
N techniques and numerical simulations are used to study various properties of the model. The scalar and 
pseudoscalar spectrum, the approach to the continuum and chiral limits, the size of the 1/iV corrections, and the 
effects of the zero momentum fermionic modes on finite lattices are studied. Also, some interesting observations 
are made by viewing the model as an embedding theory of the Higgs sector. 



1. Introduction 

When the high frequency modes of the gauge 
and fermionic fields of QCD are integrated down 
to the energy scale E corresponding the correla- 
tion length of the gauge field (E « glueball mass 
1550 MeV (l]) the resulting effective theory 
will essentially be a theory of fermions with con- 
tact interactions and cutoff A < 1550 MeV. The 
resulting effective Lagrangian will maintain the 
original chiral symmetry but will be more compli- 
cated. At energies much smaller than A it should 
be enough to keep in the Lagrangian the least ir- 
relevant operator, namely the four-Fermi dimen- 
sion six operator. This is one way [^| to motivate 
the Nambu-Jona-Lasinio (NJL) model. 

Unfortunately, by only keeping the four-Fermi 
operator, valuable information is lost and the 
NJL model does not confine the quarks. Further- 
more, if, for example, we want to study the a 
particle, which on phenomenological grounds is 
believed to have mass ~ 750 MeV, then the sep- 
aration of scales is probably not large enough to 
justify the neglect of operators with dimension 
higher than six. Nevertheless, the NJL model 
possesses the same chiral symmetry as QCD and 
it can realize this symmetry in the Goldstone 
mode. It is this feature that is most crucial in 
the understanding of the lightest hadrons. Fur- 
thermore, our interest in the model is not so much 
aimed at its quantitative predictability but rather 
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on the qualitative insights that can provide as a 
low energy theory of QCD, as an embedding the- 
ory of the Scalar Sector of the Minimal Standard 
Model and as a four-dimensional interacting the- 
ory of fermions on the lattice. 

The NJL model has been studied extensively 
for various cases with continuum type regulariza- 
tions. For a comprehensive review the reader is 
referred to || and references therein. Further- 
more, the NJL model is a special case of Yukawa 
models that, under a different context, have been 
studied extensively with lattice regularization [|J . 
The model has also been studied on the lattice 
|i| in connection with the possible equivalence of 
the top quark condensate with the Higgs field || . 
Also in this conference work was presented re- 
garding the study of the NJL model on the lattice 
with staggered fermions in connection with QED. 

In this work Q we consider the two flavor (up 
and down) NJL model with SU(2) x SU(2) chi- 
ral symmetry and SU(N) color symmetry, with 
scalar and pseudoscalar couplings || on the four- 
dimensional hypercubic lattice and we consider 
both naive and Wilson fermions. We study the 
NJL model using large N techniques and obtain 
analytical results both on finite and infinite vol- 
umes. The infinite volume results are obtained 
sufficiently close to the continuum limit using 
asymptotic expansions. We also study the model 
for N = 2 using an HMC numerical simulation 
ficl l with Conjugate Gradient and leap-frog algo- 
rithms. 
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2. Results 

The seven main results that stem from our 
analysis are presented below. 

1) For naive fermions we calculate at large N 
and with M„ = 140 MeV the a mass (M CT ), the a 
width (Tcr) and the constituent quark mass (M q ) 
in physical units as functions of the cutoff. By 
setting M q = 310 MeV, we find M a = 726 MeV, 
r CT = 135 MeV, and A = n/a = 1150 MeV. M a 
is consistent with phenomenological expectations 
and A is consistent with the expectation that the 
cutoff should be close and below the mass of the 
lightest glueball (1550 MeV). The width however 
is underestimated. The reason is traced to the 
fact that to leading order in large N the width 
receives contributions only from the quark bub- 
ble and not from the pion bubble because the pion 
bubble is of order \/N . Because the phase space 
available for the a to decay to two quarks is much 
smaller than the phase space to decay to two pi- 
ons the pion loop contribution, although of order 
l/N, is probably more important than the quark 
loop contribution. 

2) The above result is relevant not only for the 
low energy QCD but also for the Higgs sector. 
It is well known that there is an equivalence be- 
tween the a-ir sector of QCD with the scalar sec- 
tor of the Minimal Standard Model. In the for- 
mer the scale is set by the pion decay constant 
(Ftt — 93 MeV) and in the latter by the weak 
scale (F^ = 246 GeV). As mentioned above we 
find that in accordance with phenomenological 
expectations in the a-ir QCD sector M^/F^ « 8, 
but in the Higgs sector all previous analysis pre- 
dicts a triviality bound of the Higgs mass with 
Mv/Ft, < 2.8 (see for example fill). In the past 
this has been a reason for concern since it could 
imply that the Higgs mass bound may be under- 
estimated. Our analysis suggests that this appar- 
ent discrepancy appears because the Higgs mass 
bounds are traditionally obtained for m a < 0.5 
while the Ma/F^ w 8 ratio is obtained for m a ks 2 
and it should therefore be accompanied by large 
deviations from the low energy behavior. Never- 
theless, this is only a suggestion since we have not 
calculated deviations from the low energy behav- 
ior of a physical process that would enable us to 



make exact statements. However, the value of the 
width serves as an indication of the size of such 
deviations. In a way, departure from low energy 
behavior will be signaled by an increasing width 
of the a to two quark decay. At m a « 2 the width 
is already fairly large. 

3) If the Higgs sector is the low energy effective 
field theory of a NJL model (which in turn is an 
effective theory of a high energy QCD-like the- 
ory) with N c = 3, rif — 2 and M„ = 0, then if we 
set the fermion mass to -p 2 - w as is the case 
for the low energy sector of QCD, the Higgs mass 
will be M a — 1915 GeV. This corresponds to 
m a = 2 and at this point one would expect very 
large deviations from the low energy behavior of 
scattering cross sections. This suggests that as 
the CM energy is turned up, first the deviations 
from the low energy behavior will become large 
signaling the onset of new physics, and later on 
the Higgs particle would be observed. 

4) With Wilson fermions we obtain at large N 
analytical expressions of the pion mass (m w ) and 
constituent quark mass (m q ) in lattice units as 
functions of the bare parameters of the model. 
We are then able to make exact statements re- 
garding the approach to the continuum and chi- 
ral limits. We draw the "phase diagram" and 
identify the single point where a continuum chi- 
ral limit (m,r — ► 0, m q — > 0) can be achieved. 
This may provide an insight on how the retrieval 
of the continuum chiral limit is achieved in QCD. 

5) At large N and for Wilson fermions the a 
particle has mass proportional to the cutoff. Our 
analysis traces this to two related reasons. First, 
although the Wilson term has raised the masses 
of the doublers to the cutoff, it has not decoupled 
them from the theory. Through vacuum polar- 
ization these contribute to the a self energy and 
raise its mass. Second, although the Wilson term 
has not altered the low frequency behavior of the 
propagating quark, it has, however, altered its 
high frequency behavior. Again, through vacuum 
polarization the high frequency modes contribute 
to the a self energy and also raise its mass. Such a 
phenomenon may also be responsible for the diffi- 
culty in observing a a particle in numerical simu- 
lations of QCD with dynamical Wilson fermions. 
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6) The numerical simulation is performed on fi- 
nite lattices. For naive fermions one would expect 
to be able to see some indication of the chiral 
phase transition as well as a a particle. However, 
by simply looking at the graph of the vacuum ex- 
pectation value vs. the coupling on an 8 3 x 16 
lattice one can not see an indication of a phase 
transition. Also the a particle in a 16 4 lattice is 
either non existent or too heavy to be measured. 
Both of these unexpected results can be explained 
at large N. The reason is traced to the existence 
of zero quark momentum modes that on a finite 
lattice are not sufficiently suppressed. The zero 
modes besides obscuring some of the physics are 
also probably partially responsible for the large 
inversion times in the HMC algorithm. To lead- 
ing order at large N the smallest eigenvalue of the 
matrix that has to be inverted is m 2 and corre- 
sponds to the zero quark momentum mode. For 
small m q , the condition number of the matrix is 

for r = and 64r/m^ for r = 1 and it is 
clear that it depends strongly on the presence of 
the zero modes. A large condition number will 
make the inversion of M'M very slow. Further- 
more the spacing of the smallest eigenvalues be- 
haves like 1/L 2 and for larger lattices the inver- 
sion times will rapidly get worse. An important 
observation can be made by noticing the depen- 
dence of the condition number on r. This suggests 
that performing the simulation with smaller r will 
yield a quite faster inversion. It is possible that 
this may also be the case for QCD. 

7) The observables measured in the numerical 
simulation (chiral condensate, vacuum expecta- 
tion value, pion wave function renormalization 
constant, and pion mass) have values that are in 
good agreement with leading order large N. This 
provides a quantitative prediction for the size of 
the l/N corrections. In agreement with the large 
N predictions discussed in (5) and (6) above, 
the a mass was very heavy to give a good sig- 
nal and was not measured. Also, measurements 
of the sigma width were not performed, but, as 
discussed in (1) above, we expect the l/N cor- 
rections to the width to be large. 

There are some interesting issues relevant to 
lattice work that have not been considered in this 
paper. It would be important to calculate the 



three and four point vertices and therefore be able 
to calculate scattering amplitudes and their de- 
parture from low energy behavior as well as the 
l/N corrections to the width. It would also be 
interesting to study the NJL model at finite tem- 
perature. Finally, it would be important to in- 
clude vector meson couplings (see, for example, 
U , |l^] ) and confirm that for the case of Wilson 
fermions the vector meson masses scale appropri- 
ately and do not become of the order cutoff as the 
a particle does. 
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